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Power corrections and renormalon resummation for the average thrust* 

Einan Gardi and Georges Grunberg a 

a Centre de Physique Theorique de l'Ecole Polytechnique 91128 Palaiseau Cedex, France 
email: gardiOcpht . polytechnique . f r , grunberg@cpht . polytechnique . f r 

Infrared power corrections for the average thrust (T) in e + e~ annihilation are analyzed in the framework of 
renormalon resummation, motivated by analogy with the skeleton expansion in QED and the BLM approach. 
Performing the "massive gluon" renormalon integral a renormalization scheme invariant result is obtained. We find 
that a major part of the discrepancy between the known next-to-leading order (NLO) calculation and experiment 
can be explained by resummation of higher order perturbative terms. This fact does not preclude the infrared 
finite coupling interpretation with a substantial 1/Q power term. Fitting the regularized perturbative sum plus 
a 1/Q term to experimental data yields af* (M z ) = 0.110 ± 0.006. 
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1. Introduction 

Power corrections to event shape observables 
in e + e~ annihilation have been an active field 
of research in the recent years. Event shapes, 
as opposed to other inclusive observables, do not 
have an operator product expansion, so there is 
no established field theoretic framework to anal- 
yse them beyond the perturbative level. On the 
other hand, the experimental data now available 
cover a wide range of scales and thus could pro- 
vide an opportunity to test QCD and extract a 
precise value of a s . 

The state of the art in perturbative calcula- 
tions of average event shape variables is 0(a 2 ), 
i.e. NLO. It turns out that experimental data are 
not well described by these perturbative expres- 
sions, unless explicit power corrections, that may 
be associated with hadronization, are introduced. 
Renormalon phenomenology allows to predict the 
form of the power terms while their magnitude is 
determined by experimental fits. 

In the work reported here |IJ we assume the 
existence of an Abelian like "dressed skeleton ex- 
pansion" in QCD and calculate the single dressed 
gluon contribution using the dispersive approach. 
This way we perform at once all order resumma- 
tion of perturbative terms which are related to the 
running coupling (renormalons) and parametriza- 



tion of power corrections. We discuss the ambigu- 
ity between the perturbative sum and the power 
corrections and show that the resummation is es- 
sential in order to extract the correct value of a s 
from experimental data. 

2. Average thrust in perturbation theory 

To demonstrate the method proposed we con- 
centrate on one specific observable, the average 
thrust, defined by 



T 



(1) 



*Talk given by E. Gardi at the QCD'99 conference, Mont- 
pellier, July 1999. 



where i runs over all the particles in the final 
state, pi are the 3-momenta of the particles and 
Ht is the thrust axis which is set such that T is 
maximized. It is useful to define t = 1 — T which 
vanishes in the 2-jet limit. 

Being collinear and infrared safe, the average 
thrust can be calculated in perturbative QCD to 
yield 

(t)^o(Q 2 ) = % [taa^(Q 2 ) + ha^(Q 2 )] (2) 

where a — a s /n : Cf = § and the coefficients 
are to = 1-58 and t\ = 23.7 — 1.69 Nf (see refs. 
inJJ). Using the world average value of a s , 
af s (M z ) = 0.117, the NLO perturbative result 
(|^) turns out to be quite far from experimental 
data. This is shown in fig. |[ From the figure it is 
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Figure 1. (t) as a function of Q = y/s for 
af* (M z ) = 0.117. Naive MS results (LO, NLO) 
and resummation results vs experimental data. 



clear that the NLO correction is quite significant, 
so the perturbative series truncated at this order 
is not very reliable. Furthermore, the renormal- 
ization scale dependence is significant. It follows 
that higher order corrections that are related to 
the running coupling cannot be ignored. 

3. Renormalons and power corrections 

In order to take running coupling effects into 
account it is useful to assume, in analogy with 
the Abelian theory, that there exists a "dressed 
skeleton expansion" . Then, the most important 
corrections which correspond to a single dressed 
gluon can be written in the form of a renormalon 
integral 

°o= J ^r^k 2 IQ 2 )a{k 2 ) (3) 

where a represents a specific "skeleton effective 
charge", not yet determined in QCD (in (l) sev- 
eral schemes were used). As opposed to the stan- 
dard perturbative approch (Q), by performing the 
integral (||) over all scales one avoids completely 
renormalization scale dependence. 

The integral (||) represents a non Borel 
summable power series. Indeed, it involves inte- 
gration over the coupling constant in the infrared 
which is ill-defined in perturbation theory due to 



Landau singularities. The ambiguous integral (||) 
can be defined mathematically, e.g. as a principle 
value of the Borel sum: Z?o|pv- This, however, 
does not solve the physical problem of infrared 
renormalons: information on large distances is re- 
quired to fix the ambiguity. The perturbative cal- 
culation contains some information about the am- 
biguity: if the leading term in the small momen- 
tum expansion of <j> is 4>{k 2 /Q 2 ) — C n {k 2 / Q 2 ) n , 
the leading infrared renormalon is located at n//?o 
in the Borel plane, and a power correction of 
the form 1/Q 2n is expected. Having no way 
to handle the problem on the non-perturbative 
level, it is natural to attempt a fit of the form 
A)|pv + A/Q 2 ". 

A stronger assumption || is that the "skeleton 
coupling" can be defined on the non-perturbative 
level down to the infrared. Then the integral 
should give at once the perturbative result 
plus the correct power term. Since the in- 
frared coupling is not known, it is considered as 
a non-perturbative parameter. Using the cut- 
off regularization of (||), one fits the data with 

A)U + W<9 2 ", where 

A>L= / ^4>(k 2 /Q 2 )a(k 2 ) (4) 
J It* K 

is fully under control in perturbation theory and 
the normalization of the power term 

rtf dk 2 

is a perturbatively calculable coefficient times a 
moment of the coupling in the infrared. Since the 
coupling is assumed to be universal, the magni- 
tude of power corrections can be compared be- 
tween different observables ||. 

The generalization of this approach to 
Minkowskian observables such as the thrust was 
discussed in ||,|l| . At the level of a single gluon 
emission it is based on a "gluon mass" renormalon 
integral f| 

r°° du 2 

R = / -4- M/* 2 ) Hv 2 /Q 2 ) (6) 

Jo M 

where the characteristic function T is calculated 
based on the matrix clement for the emission of 
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one massive gluon and a off is related to the time- 
like discontinuity of the coupling. Regularizations 
of the perturbative sum (61) with a at one or two 
loops were discussed in [jj. We skip it here for 
brevity. 

4. Fitting experimental data 

In order to perform renormalon rcsummation 
at the level of a single gluon emission, we calcu- 
lated the characteristic function for the thrust [Q . 
Then Rq is evaluated taking either the one or two 
loop coupling. Since R does not exhaust the full 
NLO correction of eq. (||), we add an explicit NLO 
correction, 5 NLO = m — t^) a 2 ^ , where t\ corre- 
sponds to (||) and tl is the piece included in R . 
It was shown in that the Abelian part of t\ 
almost coincides with that of t\ in spite of the 
non-inclusive nature of the thrust. This is crucial 
for the applicability of the current resummation 
approach. Finally, we add an explicit power cor- 
rection X/Q 2n = X/Q where the power n = \ 
is determined from the small fj, 2 expansion of T. 
Thus, using the PV regularization of Rq, we fit 
the data with 
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The resummation as well as the fit results are 
presented in fig. [I] for the world average value of 
a s . The resummation by itself is quite close to 
the data. Note also the stability of the result as 
a is promoted from one to two loops. Next, we fit 
also the value of a s . The best fit % 2 /point = 1.35, 
obtained with of* (M z ) = 0.110 ±0.006 and A = 
0.62 ± 0.11, is shown in fig. |. 

5. Truncated expansions in MS 

It is interesting to compare the resummation 
to a truncated expansion in a-^ ■ Expanding Rq 
with a at 1-loop we obtain a series of the form 
E^o^olt 1 . The coefficients t° can be consid- 
ered as predictions for the perturbative coeffi- 
cients fj, provided the coupling a is close to the 
correct "skeleton coupling" , and provided that 
the "leading skeleton" Rq is indeed dominant al- 
ready in the sub-asymptotic regime. The signifi- 
cance of the NNLO and further sub-leading cor- 



rections which correspond to the dissociation of 
the emitted gluon is demonstrated in fig. g. The 
figure shows 
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where k is the order of truncation for k = 2 
through 6. At these orders the series is still con- 
vergent. We conclude that a major part of the 




Figure 2. (*) per t as a function of Q for 
a" s (Mz) = 0.111. The perturbative series is 
gradually improved by adding terms from Rq. 



discrepancy between the NLO result and experi- 
ment is due to neglecting this particular class of 
higher order corrections. Next consider a fit to 
experimental data based on the truncated expan- 
sion of Rq (§): (t) + X peT t/Q. The fit results 
are listed in table EL 

Table 1 

Fit results based on the truncated expan- 
sion in MS 



k 


af s (M z ) 


Apcrt (GeV) 


2 


0.128 


0.72 


3 


0.118 


0.65 


4 


0.115 


0.58 


5 


0.114 


0.50 


6 


0.114 


0.40 


PV 


0.111 


0.73 
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The quality of the fit is roughly the same in all 
cases. We see that the resummation is absolutely 
necessary in order to extract a reliable value of a a . 

6. Infrared cutoff regularization 

Putting an infrared cutoff /ij on the space-like 
momentum (||) we separate at once the perturba- 
tive and non-perturbative regimes as well as large 
and small distances. The cutoff regularized sum, 
Ro\in f or P"i = 2GeV is shown in fig. || together 
with the PV regularization. For the lower Q val- 



2— loop resummation + NLO 

PV regularization 

cutoff regularization 

+ VQ fit 




Figure 3. (t) as a function of Q for 
af s (Mz) = O.f fO. Resummation in different reg- 
ularizations and a fitted curve based on either of 
them. 



ues the difference between the two regularizations 
is large. This means that a large contribution to 
i?o|p V comes from momentum scales below 2 GeV, 
where the coupling is not controlled by perturba- 
tion theory. Finally fitting the data 



(*) 



C f 



x 



ILL 

Q 



(9) 



we get the same result as with the PV regular- 
ization (0) . The only difference is in the required 
power correction A P/ = 1.49. This is general: 
the two regularizations differ just by (calculable) 
infrared power corrections - in our case 1/Q - 
and are therefore equivalent once the appropriate 
power terms are included. 



7. Conclusions 

The assumption of a "skeleton expansion" im- 
plies that resummation of perturbation theory 
and parametrization of power corrections must 
be performed together. For the thrust renor- 
malon resummation is significant and closes most 
of the gap between the standard perturbative re- 
sult (NLO) and experiment. The resummation is 
crucial to extract the correct value of a s . 

The infrared sensitivity of the thrust leads to 
ambiguity in the resummation, which is settled 
by fitting a 1/Q term. In the infrared cutoff reg- 
ularization this power term is substantial. 

The resummation leads to a renormalization 
scale invariant result. In the BLM approach it 
corresponds to a low renormalization point in MS, 
MIL- 0.0447 Q. 

E. De Rafael What is the physical meaning, 
in QCD, of the scale of the 1/Q correction? What 
other processes are sensitive to these corrections? 

The first question is basically an open one. 
Deeper understanding could hopefully be gained 
once renormalon phenomenology is supported by 
more rigorous field theoretic methods. In the in- 
frared finite coupling approach the 1 /Q correction 
is understood as a moment (j^) of a universal in- 
frared finite coupling. This allows comparison be- 
tween observables - for example the C parame- 
ter is sensitive to similar 1/Q corrections. In the 
framework of shape-functions J^/, a relation with 
the energy-momentum tensor was suggested. 
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